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Vortex configurations in rotating Bose-Einstein condensed gases trapped in power-law and an- 
harmonic potentials are studied. When the confining potential is steeper than harmonic in the 
plane perpendicular to the axis of rotation, vortices with quantum numbers larger than one are 
energetically favorable if the interaction is weak enough. Features of the wave function for small 
and intermediate rotation frequencies are investigated numerically. 



I. INTRODUCTION 

' Trapped Bose-Einstein condensed gases provide a novel kind of condensed matter system which can sustain quan- 
tized vortices, as has been realized in the recent years ^, ^ Features which differ from those in macroscopic 
superfLuid systems are expected here. The fact that these systems are finite, and that the parameter regime where 
the vortex core size is comparable to the system size is attainable, in principle opens up for the possibility for vortices 
in these systems to have a circulation quantum number q larger than unity. 

It is well known that multiply quantized vortices are not thcrmodynamically stable in spinless, macroscopic and 
• homogeneous superfluids, because the energy of a vortex depends on the square of the circulation, and therefore two 
singly quantized vortices with a finite spatial separation have lower energy than one doubly quantized, while giving 
04 1 the system the same angular momentum However, multiply quantized vortices may well be energetically favorable 
' in systems which do not fulfill the criteria of being homogeneous, spinless and macroscopic. The assumption of 
homogeneity is not met in superconductors with pinning forces, where indeed doubly quantized vortices are observed 
In ^He-A, where the order parameter is not a scalar due to spin degrees of freedom, a lattice of doubly quantized 
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vortices with filled cores has recently been observed Q]. The argument also fails in spatially confined systems where 
■ the vortex cores are not much smaller than the system, such as mesoscopic superconducting disks, where vortices with 
large quantum numbers have been predicted |Q . 

In condensed Bose gases confined in harmonic-oscillator potentials, it has been found both analytically numer- 
ically Ip^ , [III and variationally that multiply quantized vortices are not energetically favorable, and no quantum 
numbers larger than unity have been observed in experiment ^. However, none of these studies has considered 
whether altering the power law of the confining potential can open up for the existence of multiply quantized vortices. 
In this paper, we shall give rigorous criteria for the thermodynamic stability of multiply quantized vortices in trapped 
Bose-Einstein condensates, and show the decisive role played by the shape of the potential. 
Q ' The system under study is a gas of bosons of mass m in an external potential T^(r'), dilute enough and at sufficiently 
O low temperature that it is well described by the Gross-Pitaevskii equation Q 
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^^m^ = -:^VV(r-, t) + V{t)^{t, t) + t/olV'lr, t)PV'(r, t) ~ nL^r, t). (1) 



' The so-called condensate wave function '0(r, t) is normalized to the number of particles N in three dimensions, and 
to the number of particles v per unit length in two dimensions. The coefficient t/o in front of the nonlinear term is 
the interaction strength, defined as Uq = Airfi^a/m, where a is the s-wave scattering length. It turns out that the 
effective measure of interaction strength in two dimensions is the product 47ri/a, which we shall denote by g. In three 
dimensions the corresponding quantity depends on the trapping potential. 

A centrifugal term is present in the Gross-Pitaevskii equation, corresponding to a rotation of the trap with the 
frequency about the z axis. At certain critical values Jlcg of the rotation frequency, we expect there to be a discrete 
transition between states of different circulation numbers, so that when ficg+i > ^ > ^cq, a state with total circulation 
quantum number q is the ground state The critical frequencies are functions of the coupling strength Uq |15| . 

When the Gross-Pitaevskii equation is valid, the total energy of the system is given by the mean-field energy 
functional 



E = J drrir) [~^^ + Vir) ~ ^Lzj i^{r) + ^C^ol^Wr- (2) 

We shall in the following be concerned with the lowest-energy solutions of the Gross-Pitaevskii equation (|l|) for a 
range of values of the driving frequency f2, interaction strength Uq and for different external potentials V . The paper 
is organized as follows. In Section we present an analytic study of systems contained in three- and two-dimensional 



power-law and anharmonic traps. In Sec. Ill, we illustrate these findings numerically for the two-dimensional case 
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and study some features of the different states with small and intermediate circulation quantum numbers. Sec. IV 
provides a conclusion. 

II. VORTEX CONFIGURATIONS IN POWER-LAW AND ANHARMONIC TRAPS 

Consider a Bose gas trapped in a cylindrically symmetric power-law potential of power n, which in cylindrical 
coordinates is written 

Vir,cl,,z)^hLo(j) +V\\{z). (3) 



The trap frequency oj, which determines the strength of the potential, defines a trap length = y^h/muj. The 
dependence of the potential on the z coordinate is totally arbitrary as long as the problem remains separable in the 
noninteracting limit, as is the case here. We shall hereafter neglect to mention the irrelevant z-dependent term and 
its associated degrees of freedom. 

In the analysis of this section, we shall consider the weakly-interacting limit where Uq is small, because if a multiply 
quantized vortex can ever be present, it will be so when the interaction is weak. To see why, note that the vortex 
cores are large for weak interactions and vice versa Therefore, in the limit of large Uq: the vortices are much 

smaller than the system size and the radius of curvature of the density profile. The result of a homogeneous bulk 
system thus applies, namely that a vortex array is the energetically favorable configuration. Therefore, we expect 
multiply quantized vortices to show up only in the realm of small (and possibly intermediate) Uq. 

Ref. |9| showed why multiply quantized vortices are never the minimum-energy state in a system which is harmon- 
ically confined in the x-y plane. We will here extend the same analysis to the case when the confinement is steeper 
than harmonic. The argument is based on perturbative treatment of the interaction term in the Gross-Pitaevskii 
equation (|l|). We therefore first study the noninteracting gas. When Uq = 0, the Gross-Pitaevskii equation is identical 
to the one-particle Schrodinger equation, whose solutions are well known when the radial confinement is harmonic. 
The eigenstates (pn^q can be labeled by a radial and an axial quantum number rir and q (as well as a quantum number 
associated with the z direction, irrelevant for our considerations). The eigenenergies including the centrifugal term 
are 

En^q = h{uj - il)q + huj{l + rir). (4) 

We shall concentrate on the states with no radial nodes, ipoq, because these have the lowest energy for a given angular 
momentum Nqh. These are in fact vortices with quantum number q. When the driving frequency fl is less than the 
trap frequency lu, the q — state is the ground state; when fl = uj, all the states ipoq are degenerate. When fl > uj, E 
decreases with increasing q, and there is no ground state. Hence, in the noninteracting case all the critical frequencies 
ricq are equal to the trap frequency uj. 

Following Ref. §, we now consider perturbatively the interaction energy. The perturbation for a given wave function 
■>p is, from Eq. (||), 

H' = luo J dr\^ir)\\ (5) 

Calculating this for general superpositions 'ip{r) = J2q^q'-Poq{f): we find that the degeneracy at f7 = w is lifted, and 
the critical frequencies flcq split apart and assume values less than uj. Fixing a total angular momentum Lz = NqU, 
it turns out that wave functions which are superpositions of states of different q have less interaction energy than the 
pure multiply-quantized vortex configurations ipoq for q > 2. Such superpositions are vortex arrays, and these are the 
energetically favorable states for any finite interaction strength. 

We shall now show how the analysis of Ref. |^ works when the trapping power n in Eq. (^) is larger than two. In 
this case, the single-particle energy varies faster than linear with the angular momentum. This implies that in the 
noninteracting case, pure multiply-quantized vortex states ipoq have lower energy than any superposition having the 
same angular momentum. Hence, it takes a finite (albeit possibly very small) interaction strength Uq to overcome this 
difference. We stress that in order to arrive at this result, we did not need to know any details of the eigenfunctions 
of a general power-law potential: the fact that the energy depends stronger than linearly on the angular momentum 
is enough. For a finite but weak enough interaction, the multiply-quantized vortex states are the minimum-energy 
configurations within their respective ranges of external rotation frequency Q. When the interaction is strong enough, 
however, the larger interaction energy of the multiply-quantized vortex states will overcome the difference in kinetic 
plus trap energy, and the ground state will again be a vortex array. This makes sense, since as we noted above, in the 
limit of strong interaction the results of the infinite and homogeneous case apply. 
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For potentials weaker than harmonic, n < 2, no rotation is possible. The potential in a frame rotating with angular 
frequency fl is V{r) — ^f^f^, which is not confining if > and V{r) is weaker than harmonic. Therefore a rotating 
state in such a trap can be at most metastable. 

If a pure power-law potential which is steeper than harmonic may seem an artificial construction which is hard to 
fabricate in practice, we note that the analysis can be repeated for an anharmonic potential, 

y(r,0,z) = imc.V(l + Aj) + V^|(z). (6) 

The important point is that for any positive anharmonicity A, the single-particle energy varies faster than linear with 
the angular momentum. This again leads to the conclusion that multiply quantized vortices can exist in anharmonic 
traps as long as A > and the coupling is small. 

We have shown that multiply quantized vortices can be stable only in trapping potentials which are steeper than 
harmonic in the x-y plane, and only for weak enough interactions. A close inspection of the above argument reveals 
the physical mechanism behind the transition. The kinetic energy gives rise to a repulsion between the vortices 
because of the overlapping velocity fields, while the trapping potential strives to keep the vortices together in order 
to reduce the spatial extent of the cloud. The interaction energy favors states with several separated vortices because 
such states have lower density on the average. The thermodynamic stability of multiply quantized vortices is thus a 
consequence of the competition between these three energies, and the power 2 is the limiting value. 



III. NUMERICAL RESULTS FOR TWO-DIMENSIONAL SYSTEMS 

We now perform a numerical study to quantify the predictions of Section ||. We shall treat the two-dimensional case 
only. We insert the power-law potential, Eq. (||), into the Gross-Pitaevskii equation (|]) and scale out the dimensions 
by defining r = dtf, t — uj^^i, ip = ^/vd^^tp, = hLz, ^ = wil, g ~ Airi/a, and for the total mean-field energy 
E = huoE. There results a dimensionless Gross-Pitaevskii equation, 

- --V^V + -r'V + fflV'l V - ^L,^, (7) 
at 2 2 

which depends on three parameters, namely the coupling g, the rotation frequency fi, and the trapping power n. We 
solve it in two dimensions using the split Fourier method fTll . By propagating Eq. in complex time, the system 
relaxes to a local energy minimum which depends on the initial state. In order to find the global minimum-energy 
configuration for each choice of parameters g, n, n, we perform the relaxation a few times using different, irregular 
initial states, compute their energy with the aid of Eq. (H), and out of the different final states we choose that which 
has the lowest energy. The grid size is 64 x 64 points, and the lattice constant is equal to 0.2 times the trap length 
dt- We have confirmed that doubling the grid size and halving the lattice constant does not alter the results. 

Figure [I] is a phase diagram for the quartic potential {n = 4), which shows how the fl-g plane is divided into 
regions of different vortex configurations with the total circulation q ranging from I up to 4. The different vortex 
configurations are labelled by the total number of quanta followed by the number of singularities present in the system: 
the configuration 3-3 thus denotes three q — I vortices whereas 3-1 is a state with one triply quantized vortex present. 
The region marked by is the nonrotating state. Since the numerical computation is done with limited precision, we 
need an operational criterion for the existence of a multiply quantized vortex. A g-fold quantized vortex is defined 
as a configuration in which the calculated circulation around a circle enclosing the origin is equal to q and only one 
density minimum can be seen on optical inspection of the plotted wave function, i. e. the inter-vortcx distance is less 
than the lattice constant of the numerical grid. The transition frequencies found using this criterion coincides with 
the point where the angular momentum saturates at an integer value. 

As expected from the analysis of the previous section, we find multiply quantized vortices for small values of g, 
which split into arrays of singly quantized vortices for larger g. The transition from a multiply quantized vortex to 
an array is continuous, as illustrated in Figure |[ where two vortices are seen to merge as the critical line is crossed. 

It is convenient to define a new quantity to describe the regime of stability of multiply quantized vortices. The g-th 
critical coupling gcq is defined so that when g > gcq, there can exist no ij'-fold quantized vortices for any value of il, 
but when g < g^q, a g-fold quantized vortex is the energy minimum within certain limits of SI. The critical coupling 
is a function of the trap power n. Fig. ^ suggests that the critical coupling increases with q. The physical reason for 
this is, that if the system is to accommodate many singly quantized vortices, the cores must be small and g large. 
Our analysis in the previous section showed that gcq = for the harmonic trap, n = 2, and gcq > when n > 2. This 
is confirmed in Figure ^, which shows how gc2 depends on n. The data points suggest that the curve gc2{n) goes to 
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FIG. 1; Phase diagram for a quartically confined two-dimensional Bose gas subject to a force rotating with an angular frequency 
Q. The state marked by is the nonrotating state, and 1-1 is a state with one central q ~ 1 vortex. The regions 2-2, 3-3 and 
4-4 denote vortex lattices containing two, three and four singly quantized vortices, respectively, while the states 2-1, 3-1 and 
4-1 are states with several quanta of circulation but only one phase singularity. 



zero linearly when n 2 from above. 

We have pursued the analogous anharmonic case, Eq. (|^), using the same numerical procedure. The qualitative 
features of the phase diagram are not different from the pure power-law case. The critical coupling gcq is now a 
function of the anharmonicity A. The result for gc2 is shown in Figure ^. The limit A — > cx) is the quartic potential, 
whose critical coupling gc2 = 30.5 can be read off from Fig. When A ^ 0, on the other hand, the potential 
approaches a harmonic one, and gc2 approaches zero. In the experimentally realized two-dimensional condensates of 
Ref. g varies between 100 and 10'*, the former for a number of particles N = 10^. Diminishing the number of 
particles by a factor 10 and introducing an anharmonicity larger than 0.1 would take these systems below gc2- 

The phase diagram in g-fl space becomes richer for higher quantum numbers q. For q < 7, there are only two 
phases, a vortex-array phase and a multiply-quantized vortex phase. For higher quantum numbers, however, a new 
kind of pattern appears, where one multiply quantized vortex is surrounded by singly quantized vortices. Figure |^ 
shows the different vortex phases for 8 < g < 11. The power of the potential is chosen to n — 2.2. We find a number 
of distinct phases, which we label by their total circulation followed by the number of singularities just as in Fig. [|. 
There are three phases containing eight quanta of circulation: a central q = 8 vortex for small g (denoted 8-1); an 
array of eight vortices for large g (8-8); and for intermediate g, a central q = 2 vortex surrounded by six q — 1 vortices 
(8-7). The transition from 8-1 to 8-7 is continuous: the six q = 1 vortices separate from the central one in a smooth 
way as the critical line is crossed. The two phases 8-7 and 8-8, on the other hand, have different symmetries and the 
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FIG. 2: Numerically computed condensate wave functions, illustrating how the distance between two vortices continuously 
decrease as the external rotation frequency il increases. The bright shades indicate high density and vice versa. The confinement 
is quartic, the coupling constant g = 28.0, and Q/uj — 1.59 for the top left panel, 1.72 (top right), 1.83 (bottom left) and 1.86 
(bottom right). The latter is the splitting frequency, operationally defined as the frequency where no mode of plotting the wave 
function reveals more than one density minimum. 

transition between these two is discontinuous. For the case of nine quanta, there exist the corresponding phases 9-1 
and 9-9, and one intermediate phase 9-8 with a central q = 2 vortex and seven q = 1 vortices. A 9-7 configuration 
with a g = 3 vortex surrounded by six q — I vortices seems to be stable in a very narrow region of width Ag ^0.2, 
between 9-1 and 9-8, but with the present numerical precision we have not been able to accurately determine the 
boundaries of this phase and it is left out of Fig. ^. For the case q=10, on the other hand, one can clearly distinguish 
four different configurations, and likewise for q = 11. 

Analogously to the q = 8 case, the transition between 9-1 and 9-7 is continuous, and so is the transition between 
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FIG. 3: Critical coupling for the stability of doubly quantized vortices, gc2, as a function of the power n for a condensate 
in a power-law trap. The analysis of Sec. ^ predicts that the curve goes to zero when the trapping power approaches 2, in 
accordance with these numerical results. 



10-1 and 10-8 and between 11-1 and 11-8. On the other hand, the configurations 9-7, 10-8 and 11-8 have entirely 
different symmetries from 9-8, 10-9 and 11-9, respectively, and those transitions are discontinuous. Finally, the three 
last-mentioned states transform smoothly into arrays of singly quantized vortices (9-9, 10-10 and 11-11) when g 
increases, by splitting of the central vortex into singly quantized ones. Figure illustrates this. 



IV. CONCLUSIONS 



We have determined the criteria for existence of vortices with circulation quantum numbers larger than unity in 
rotated trapped Bose-Einstein condensates. We have proven that in both two- and three-dimensional systems, multiply 
quantized vortices are the energy minimum if the trapping potential in the plane perpendicular to the axis of rotation 
is steeper than harmonic, and the interaction is sufficiently weak so that the vortex cores are not much smaller than 
the cloud, and the external rotation frequency is within appropriate limits. For stronger interactions, the multiply 
quantized vortices break up into arrays of several vortices. The physics of the transition is the interplay between 
the inter-vortex repulsion and the trapping potential. In the case of two-dimensional systems, we have numerically 
determined the regimes of thermodynamic stability of different vortex configurations. For clouds rotated at large 
angular velocities, there are regimes where arrays containing both singly and multiply quantized vortices are stable. 

Most trapping potentials used in experiments are harmonic, but using Laguerre-Gaussian beams one can construct 
optical traps of arbitrary even power laws, where our predictions can be tested p9[. A two-dimensional Sodium 
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FIG. 4: Critical coupling gc2 for a condensate trapped in an anharmonic potential as a function of the coefficient A of the 
quartic term. When A approaches zero, the trap approaches the harmonic case, and the critical coupling drops to zero. 



condensate can host doubly quantized vortices if the number of particles is of the order 1000 or smaller, and a trap 
anharmonicity of about 10% is introduced. Alternatively, the appropriate weak-coupling regime can be attained 
by using Feshbach resonances to diminish the scattering length In a three-dimensional system, a quantitative 
determinition of the conditions has been prohibited by numerical limitations. The technique of vortex imaging by 
interference [pit , which directly images the phase of the condensate, promises to be a convenient way to detect multiply 
quantized vortices. Alternatively, one can look for multiply quantized vortices by measuring the angular momentum 
of the system |23], or simply by looking for anomalously large vortex cores. 
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